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There are strong correlations between the three structural properties of elliptical galax- 



h 

ies - stellar mass, velocity dispersion and size - in the form of a tight "fundamental 
plane" and a "scaling relation" between each pair. Major mergers of disk galaxies are 
assumed to be a mechanism for producing ellipticals, but semi-analytic galaxy for- 
fNj ■ mation models (SAM) have encountered apparent difficulties in reproducing the 

J> ' observed slope and scatter of the size-mass relation. We study the scaling relations 

, of merger remnants using progenitor properties from two SAMs. We apply a simple 

merger model that includes gas dissipation and star formation based on theoretical 
considerations and simulations. Combining the SAMs and the merger model 
allows calculation of the structural properties of the remnants of major 
mergers that enter the population of elliptical galaxies at a given redshift. 
Without tuning the merger model parameters for each SAM, the results roughly match 
the slope and scatter in the observed scaling relations and their evolution in the red- 
shift range z — — 3. Within this model, the observed scaling relations, including the 
tilt of the fundamental plane relative to the virial plane, result primarily from the 
decrease of gas fraction with increasing progenitor mass. The scatter in the size-mass 
relation of the remnants is reduced from that of the progenitors because of a 
correlation between progenitor size and gas fraction at a given mass. 
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1 INTRODUCTION 

The merging of disk galaxies is one of the main hypoth- 
esized mechanisms for the formation of elliptical galaxies. 
Simulations have shown that mergers of disks with simi- 
lar masses can effectively disrupt the ordered rotation in 
the disks and convert it into random velocity support, cre- 
ating merger rem nants that appear simil a r to observed 
elliptical galaxies jToomre fc Toomrei Il97j; iToomrd I l977l : 



iBarnes fc Hernquisu Il992l : iMihos fc Hernquisj 11994 ) ■ Fur 
thermore, the ACDM cosmology predicts the hierarchical 
buildup of galaxies through a sequence of mergers. These 
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results suggest that merging is a likely mechanism for the 
production of elliptical galaxies. 

However, observed ellipticals follow a number of scaling 
relations, including relatively tight relations between stellar 
mass and velocity disp ersion, the Faber- Jackson relation 
IIFaber fc Jackson LL976), and between size and stellar 
mass ( Kormendv ll9T7l 'l . Furthermore, observed ellipticals 
fall on a tight plane, the fundamental plane (FP), in the 
three-dime nsional space of stellar mass, size, and velocit y 
dispersion (|Diorgovski fc Davislll987l ; iDressler et"al]|l987l '). 
Recent studies of the Sloan Digital Sky Survey (SDSS) 
have provided excellent statistics on these s caling relations 
in the local universe JShen et al.1 |2003|; I Bernardi et al.1 
l2003al lbl: iPadmanabhan et al.1 12004 : iGallazzi et al.1 l200q 
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IShankar et al]|2010ri , and studies using high redshift surveys 

have provided evidence for the evolution of these relations 
i it 1 | 

over cosmological time I Bardcn et al. 2005; Mcintos h et al 
2005 1; iTruiillo et alJlioOfj; iTruiillo et alj |2007l; ICimatti eta! 



2008 



Saracco et al.l 



van der Wei et al l 12008: 



2009 



van Dokkum ct al 



Here we implement a simplified version of this model us- 
ing post -processing of merge r outputs from the SAMs devel- 



2008 



Buitrago et al.l 

120101 ; iFan et all 



20101 ; iMancini et alll2010l ; I Williams et al.|[2oTof ). If mergers 
are a major mechanism for the production of elliptical 
galaxies then they must be able to produce the correct 
scaling relations as well as the evolution of these scaling 
relations over time. Theoretical high-resolution, hydro- 
dynamical studies have shown that simulations of gas-rich 
galaxy mergers are capable of reproducing the observed 
scaling relations of ellipt ical galaxies if the correct progeni- 
tor p r operties are used (iDekel fc Co x 2006; Robertson et al.l 
120061 ; iHopkins et"aH |200S| ; iBournaud et al.1 120111 ). These 
studies have been successful in creating high- 
redshift compact ellipticals from major mergers 
of gas-rich compact d isk galaxies (|Bournaud et aU 
2010; IWuvts et"aT]|2010h . This is a step toward verifying 
the production of scaling laws through mergers. However, 
current computing power only allows the simulation of rel- 
atively small numbers of mergers, and the space of possible 
merger initial conditions and progenitor properties is quite 
large. Furthermore, these simulations are not placed within 
a cosmological context, making it more difficult to explore 
in detail the origin and evolution of scaling relations. 

Currently, the primary theoretical tool for study- 
ing the evolution of statistical samples of galax- 
ies over cosm ological timescale s is s e mi-analytic mod- 



Somerville & Primacklll999l; ICole et al.1 200d; lHatton et al. 


2003 


ICroton et al. 20061: De Lucia et al. 20061: Bower et al. 


200£ 


Somerville et al.1 12008|; iFontanot et al.1 120091; Benson 


201C 


Benson fc Bowed [2010j; ICook et al.l 
. These models combine dark matter 


l20ld; iGuo et al. 


201C 


halo merger trees 



with analytic recipes for populating the halos with galax- 
ies. However, semi-analytic models (SAMs) do not currently 
incorporate realistic formulas for predicting the properties 
of the remnants of galaxy mergers including the effects 
of dissipation. The agreement between current SAMs and 
observed early-type scaling relations is not impressive. In 
particular, the observed size-mass relation of early-types is 
steeper than that of their potential late-type progenitors 
(|Shen et al.1 l2003h , and the scatter in the observed size- 
mass relation for early- types is remarkably small |Shen et al.l 
120031 : iNair et al.1 120101 ). The disspationless merger models 
currently implemented within SAMs h ave thus far been un- 
able to reprodu ce these features (e.g.. [Shankar et al.l 120101 : 
IGuo et alj|201fjh. Using a s i mple power law dissipation 
model iKhochfar fc Silkl {2006) were able to repro- 
duce the redshift-size evolution of elliptical galax- 
ies, which, along with the high-resolution hydrody- 
namical simulations described above, suggests that 
dissipative effects may play an important role in de- 
termining elliptical galaxy scaling relations. 

ICovington et al.1 (|2008T ) recently developed a physically- 
motivated analytic model for predicting the stellar half-mass 
radii and central velo city dispersions of merger remnants 
l|Covington et al1l2008l ). The parameters in this new merger 
model were calibrated using a suite of high-resolution hy- 
drodynamical galaxy merger simulations (see Section [2]). 



oped bv lCroton et al.l 



ulation (|Springel et al 



2006), base d on the Millennium sim - 
20051 ) . and ISomerville et al.l l|2008h . 



This results in a population of tens of thousands of merger 
remnants over a large range of redshifts (0 < z < 3) com- 
plete with predicted values of size, stellar mass, and veloc- 
ity dispersion. Comparison of the modeled population of 
ellipticals with the observed scaling relations provides an 
important test of the merger hypothesis as well as physi- 
cal insight into the origin and evolution of these relations 
via merging. It is important to note that the model 
presented here is calibrated only for major merg- 
ers of disk-dominated galaxies, and that we make 
no attempt to model subsequent evolution following 
the formative major merger. Future work will expand 
our model to include minor mergers and mergers 
of bulge-dominated galaxies. We will also implement 
our merger model self-consistently within the semi-analytic 
machinery rather than by post-processing. 

Section 2.1 describes our analytic merger model for 
calculating the properties of stellar spheroids from galaxy 
mergers including energy losses from dissipation. Section 2.2 
explains how we implement ou r mer ger model using output s 
from the ICroton et aU l|2006h and ISomerville et al.l (|2008h 
SAMs. Section 3 systematically explores the effects of the 
merger model. In order to turn the rather shallow size-mass 
relation of disk galaxies into the steeper size-mass relation 
of observed early-type galaxies, we find that it is essential 
to include both dissipation and the decreasing gas content 
of more massive progenitor disk galaxies. In Section 4 we 
summarize the observational results that we compare with 
our model outputs, focusing especially on Sloan Digital Sky 
Survey (SDSS) data for nearby galaxy size vs. mass and 
other scaling relatio ns, and on data from several surveys 
l|Truiillo et al.l l2006) for the evolution of these relations to 
higher redshifts. In section 5 we use the predicted properties 
of progenitor disk galaxies from the two semi-analytic mod- 
els to predict the size-mass and other relations for spheroids 
formed in major gas-rich mergers and compare them with 
the observations out to redshift z — 3. Section 6 summa- 
rizes these results and discusses their implications and some 
follow-on studies that are in progress. Finally, an Appendix 
describes an improvement in the treatme nt of the central 
dark m atter in the anal ytic merger mode l of lCovington et all 
(2008) that we used in lCovingtonl (|2008h and in the present 
paper. 



2 METHODOLOGY 

We use a combination of modeling approaches to construct a 
theoretical framework for predicting the evolution of early- 
type scaling relations over cosmological time. In previous 
work, a large suite of hy drodynam i cal galaxy merger sim - 
ulations were developed (|CoxN2004l : ICox et al.l [20061 . l2008r ). 
These simulatio ns were performed using the N-body/SPH 
code GADGET ijSpringel et al.ll200ll ; ISpringel fc Hernauistl 
2003), and include hydrodynamics, star formation, and stel- 
lar feedback. The simulation suite contains mergers with 
a wide variety of progenitor properties, mass ratios, and 
merger orbits. Variations in progenitor properties include 
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a range of stellar masses, gas fractions, dark matter halo 
concentrations, bulge fractions, baryonic fractions, and gas 
disk sizes. 

In subsequent work, Covington et al. (2008, hereafter 
C08) constructed a physically-motivated analytic galaxy 
merger model capable of predicting the half-mass radii, stel- 
lar masses, and velocity dispersions of galaxy merger rem- 
nants given the properties of the disk progenitor galaxies 
and the initial orbits of the mergers. This model was cali- 
brated using the galaxy mer ger suite described above. Un- 
like p revious similar models (|Cole et al . 2000; H atton et al.l 
I2003T I. this model includes the effects of star formation and 
energy loss due to dissipation. Here we combine a simplified 
version of this new galaxy merger model with merger rates 
and progenitor properties predicted by two semi-analytic 
models (SAMs) in order to explore the creation and evo- 
lution of the scaling relations of early-type galaxies. 

2.1 Description of the Merger Model 

Previous models for predicting the sizes and velocity disper- 
sions of galaxy merger remnants employed a c ombination 
of energy conservation and the virial theorem (jCole et al.l 
l200d ; lHatton et al.ll2003n . Assuming homology between pro- 
genitors and remnants then allows a straightforward calcu- 
lation of remnant properties. However, the most common 
mergers at higher redshifts are 'wet' mergers of gaseous 
galaxies, and the approach used previously does not ac- 
count for energy losses due to gas processes. Recent work 
has suggested that these dissipative effects play an impor- 



tant role in the formation of elliptical galaxies JCox et al 



2ooe 


Robertson et al.l 


2007 


Cox et al.ll2008l: 



20061; iDekel fe Coxll2006l ; ICiotti et al 



Hopkins et al.l l2008). though there 



is also evidence that disks may preserved in ma- 
jor, gas-rich me rgers (Springel &; Hernquist 2005; 
[Robertson et al.ll2006l : iHopkins et alJl2009h ~C08 de- 
velops a framework for modifying the energy conservation 
approach to include dissipative losses due to gas processes. 

In order to account for dissipative losses, we include a 
radiated energy term within the energy conservation equa- 
tion. Furthermore, since stars form during the merger, we 
include the mass of gas that will form stars when calcu- 
lating the internal energies of the progenitor galaxies. C08 
found that the energy lost and the number of new stars 
formed are a function of both the gas fraction and the initial 
merger orbit. Closer, more disruptive encounters will result 
in greater total energy losses and larger numbers of stars 
forming. However, the set of orbital parameters in the sim- 
ulations used to calibrate the model was not intended to be 
statistically representative of mergers in the real Universe. 
Rather, it was intended to sparsely span the range of rea- 
sonable values. In the present work, we find that when the 
merger model of C08 is applied to SAMs using the distribu- 
tion of merger orbit al parameters observed in cosmological 
N-body simulations dBenso n 2005) there is relatively small 
scatter in remnant properties as the result of orbit. Con- 
sequently, for this paper we use a simplified version of the 
model in C08 where remnant properties are not a function of 
orbit. This simplified model is calibrat ed to the same suite of 
galaxy merger simulations (|Coxll200"3 l. effectively averaging 
over the range of orbits within the suite. 

Additionally, the model for calculating remnant central 



dark matter fractions developed in C08, when applied within 
SAMs, fails for a subset of cases with large radii and high 
gas fractions. Therefore, here we use a new model for cal- 
culating central dark matter fractions that has a broader 
range of applicability (see Appendix A for details). To- 
gether, the modifications described result in a simpler model 
that captures the core dynamics of the more complicated 
model of C08. This has the added benefit of reducing the 
number of unknown model parameters. However, the scal- 
ing relations produced by the two models are not identical, 
and further work now underway implementing these models 
within SAMs will be needed to determine whether the sim- 
pler model presented here is sufficient (see Section 5.3 for a 
discussion of differences). 

The first step in the new model is to calculate 
the star formation effi ciency of the merger, e. Following 
Som erville et al.l |200l| ) we use 



ei:l 



Ms; 



M prin 



(1) 



This form u la wa s calibrated to our merger simulations in 
ICox et aD (|2008l ) resulting in parameter values of ei,i = 
0.55 and 7 = 0.69. For each merger M sa t is taken to be the 
total mass (dark matter plus baryons) of the less massive 
progenitor and M pr i mary is taken to be the total mass of the 
more massive progenitor. For each progenitor 

Af ns = eM gas , (2) 

where M ns is the mass of new stars formed in the merger 
and Mgas is the initial (cold) gas mass of the progenitor. 

As in C08, we use an energy conservation equation to 
calculate the half-mass radius of the merger remnant. This 
equation includes final and initial internal energy terms and 
a term for energy loss due to gas dissipation, 



-Efinal — -Emit + -B ra d- 

The energy terms are calculated as follows: 

( (M s ,l + A/ns.l) 2 | (Afs,2 + M ns , 2 ) 



-Emit = G 



-Efi na l = G 

and 

-Erad — Cr a d,/"g Ei n it 



V" #1 R2 

(M s ,i + M ns ,i + M s , 2 + M ns , 2 ) 2 



(3) 

(4) 
(5) 

(6) 



where M s ,i is the initial stellar mass of progenitor i, Ri is 
the initial 3D stellar half-mass radius of progenitor i, and 
Eflnal is the 3D stellar half-mass radius of the remnant. 
C ra d = 2.75 is a constant parameter calibrated using the 
merger simulations, and 

M g ,i + M g , 2 



M s ,i + M g ,i + M s , 2 + M, 



g,2 



(7) 



is the merger gas fraction where M Sl s and M gl i are the 
stellar and gas masses, respectively, of the progeni- 
tors. Unlike in C08, here orbital energy is neglected, 
because we find that, using orbital parame ters ob- 
serve d in cosmological N-body simulations (|Bensonl 
l2005h . the scatter resulting from orbital effects is 
small. In the limit of no dissipation (i.e./ g = 0), Eqs. 
3-5 in our paper reduce to Eqs. 1 and 2 of Naab 
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Figure 1. Half-mass radii of merger remnants from our simulation suite plotted against radii predicted from our simplified analytic 
merger model, (a), and the full C08 model with central dark matter density calculated as in Appendix A (c). Panels (b) and (d) show 
the predictions for velocity dispersions for the simple model and C08 plus Appendix A, respectively. Different symbols represent different 
types of progenitor galaxies, described in detail in C08. This new, simpler merger model does not capture variations due to initial 
orbit, which is seen in the horizontal spread within symbols of the same size, shape, and color, but we find that the spread 
resulting from this mechanism is insignificant when drawing orbits from the distribution observed in N-body simulations. The simple 
model was calibrated with a few more simulations that were not available during the initial calibration in C08. 



et al. (2009), if one assumes the proportionality of 
effective and gravitational radii as we do here. 

Once we have calculated the final radius, we can com- 
pute the velocity dispersion. First, we calculate the central 
dark matter fraction of the merger remnant, /d m ,f (see Ap- 
pendix A for details). Then, as in C08, we calculate the 
remnant velocity dispersion using the virial relation, 



M Sj f is the stellar mass of the remnant, and C s i g = 0.15 is a 
constant parameter fit to the simulations. 

In order to determine the parameters C ra d and C v i r , we 
calculated a best-fit of the model to the merger remnants in 
our simulation suite. The fit of the new model is depicted 
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in Figure []J,-b, while the results from the C08 model with 
the modified calculation of central dark matter fraction are 
shown in panels c and d. On average, the new model cor- 
rectly predicts the remnant sizes and velocity dispersions. 
However, there is more spread in the accuracy of predic- 
tions in the simple model. This results from variations in 
initial orbit in our simulation suite. As stated above, we 
find that this variation is insignificant when mergers 
are calculated for the distribution of orbits found in 
cosmological N-body simulations. 

We note that our merger model has primarily 
been calibrated for major mergers of disk-dominated 
galaxies. An expansion of this model that includes 
a detailed treatment of minor mergers and mergers 
of elliptical galaxies is forthcoming (Porter et al., in 
prep), but here we do not attempt to model sub- 
sequent evolution of elliptical galaxies following the 
initial major merger because the current versions of 
the SAMs do not track bulge sizes using our dissi- 
pative model. In limiting our model to major merg- 
ers, we are only examining one a venue of elliptica l 
gala xy formation. Re cent works (Na ab et alj 2006; 
Bournaud et al. I 120071 ) have suggested that ellipti- 
cals can also form as a product of minor mergers, 
but we do not attempt to model this here. 

2.2 Implementation of Merger Model with SAMs 

In order to make comparisons with observational data, and 
to explore the evolution of scaling relations over cosmologi- 
cal time scales, we apply the new merger model to progen - 
itor properties from two SAMs: the lSomerville et al.l (|2008r i 
SAM (here after S08) and a SA M based on the Millennium 
Simulation (|Croton et al.ll2006l ). We implement the merger 
model externally (i.e., by post-processing) rather than in- 
corporating it within each of the SAMs. This provides an 
expedient means of exploring the properties of new ellipti- 
cal galaxies arriving into the population via the merging of 
disk-dominated galaxies. We restrict our analysis to merg- 
ers of disk galaxies because the current generation SAMs 
do not have a reliable method for calculating bulge sizes, 
and we need stellar radii in order to calculate initial inter- 
nal energies. Furthermore, since the merger suite was 
designed to model mergers of disk galaxies, there 
were few bulge-dominated progenitors used to cal- 
ibrate the merger model. The performance of the 
model for mergers of bulge-dominated galaxies has 
not been tested. Consequently, we only analyze mergers in 
each SAM where both progenitors have a stellar disk more 
massive than their bulge. Additionally, since we are compar- 
ing the remnants to observations of elliptical galaxy popu- 
lations, we only include major mergers with a mass ratio of 
1:3 or greater, as these are the mergers expected to create 
elliptical galaxies. 

Each SAM directly provides the stellar mass, disk ra- 
dius, and mass of cold gas for each progenitor. In addition, 
the model requires information about the dark matter halo 
of each progenitor galaxy. Specifically, we need to be able 
to calculate mass as a function of radius in order to cal- 
culate both the halo half-mass radius and the central dark 
matter mass. S08 and Millennium provide the masses of the 
dark matter halos, but specify the halo mass distributions 



using different quantities. S08 provides halo concentrations. 
The Millennium SAM provides Vmn* and VW, from which 
concentration can be calculated. For both cases the concen- 
tration, virial mass, and redshift are used to calculate the 
distribution of mass within the halo. For S08 the spherical 
top-hat collapse model is used to c alculate the virial overden - 
sity, using the approximation from lBrvan fc Normanl l|l998T l. 
whereas for Millennium the virial overdensity is assumed to 
be 200. All of these calculations assume the ACDM concor- 
dance cosmology, with f2 m = 0.3, Qa = 0.7, and h — 0.7. 



3 SYSTEMATIC EXPLORATION OF THE 
MERGER MODEL 

To gain a better intuitive grasp of the behavior of the 
merger model, we systematically explore the effects of vari- 
ations in the progenitor properties and model parameters. 
For this study we introduce a series of four idealized progen- 
itor galaxy models. In order to isolate the effects of various 
progenitor properties we begin with a fiducial progenitor 
and scale the progenitor properties in such as way as to 
keep baryon fraction, average density, and halo concentra- 
tion constant. 

The fiducial ga laxy with which we begin our series is 
the G3 galaxy from lCox et alj (|2008T l . whose properties are 
designed to fit the average properties of observed nearby disk 
galaxies in the Sloan Digital Sky Survey. The G3 has a dark 
matter halo mass of 1.1 x 1O 12 M0, a stellar mass of 5.0 x 
1O 1O M0, an initial half-mass radius of 3.8 kpc, and a halo 
concentration of 6.0. Each subsequent galaxy in the series 
is created by reducing the mass by 1/3, keeping the average 
densities, concentration, and baryon fraction constant, and 
calculating a using the virial theorem. 



3.1 Effect of Varying Gas Fraction with Mass 

Previous work has shown that gas can have a signif- 
icant e ffect on merger r emna nts (Barnes fc Hernquistl 
19961: iRobertson et all 120061: IPekel fc Cox! 120061: 



Springel fc Hernquisdl2005l ; iNaab et al.ll2006l '). lDekel fc Cox! 
(2006) suggest that a systematic variation in gas fraction 
with mass could be responsible for th e tilt in the fun d amen - 
tal plane of elliptical galaxies, and Hopk et all (|2008h 
demonstrate observational evidence for this hypothesis. 
Thus, we examine the effect of progenitor gas fraction, 
defined here as the ratio of gas to stellar mass, on the 
properties of equal-mass merger remnants. 

In our first experiment we set the gas mass of each pro- 
genitor such that gas fraction is constant as a function of 
mass. We run the model for several gas fractions, setting 
the ratio of gas mass to stellar mass to 0.0, 0.25, 0.5, 1.0, 
and 2.0. The size-stellar mass relation, stellar-mass Faber- 
Jackson relation (F J) , and virial projection of the fundamen- 
tal plane (FP) of the merger remnants are shown in Figure 



2(a) 



Several observations are worth noting about the size- 
mass relations plotted in Figure [2(a) | First, for all gas frac- 
tions the remnant relations (crosses) are just shifted hori- 
zontally and vertically from the progenitor relation (stars) 
without any significant rotation (i.e., change of slope). That 
is, the vector between progenitor and remnant is constant 
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(a) (b) 

Figure 2. Scaling relations for 1:1 merger remnants produced by our model using a series of idealized progenitors of different 
mass. The progenitors used in (a) have a range of initial gas fractions, which are held constant as a function of mass. For remnants, 
colors depict ratio of gas mass to stellar mass in the progenitors that created that remnant where black=0.0, light blue=0.25, green=0.5, 
blue=1.0, and red=2.0. Dotted lines connect remnants resulting from progenitors with different mass but the same gas fraction. The 
progenitors used in (b) have gas fractions that vary as a power law of baryonic mass (M ga s/Mbaryons <x -^baryons ) • ^ or remnants, colors 
depict the value of 7 used to set the progenitor gas fractions black=0.0, red=0.5, and blue=1.0. Dotted lines connect remnants with 
progenitors from the same gas fraction power law. The top panels show the size-stellar mass relation of the merger remnants (crosses) 
and progenitors (stars), while the middle and lower panels show the stellar mass FJ relation and the fundamental plane for the merger 
remnants, respectively. The solid line plotted in the lowest panels shows the slope for galaxies with a virial scaling, whereas the dashed 
line shows the observed tilt in the FP. Steepening in the size-mass relationship and tilt in the FP are produced by a decreasing gas 
fraction with increasing baryonic mass. 



for any given gas fraction. This means that Rt/Ri and 
A/star.f /-Mstar,j ar e also constant for any particular gas frac- 
tion. This is worth noting because the observed size-mass 
relations for disks and ellipticals are significantly rotated 
from one another with the relation for ellipticals being much 



steeper. A realistic mechanism for the production of ellipti- 
cals must account for this steepening. 

Remember that the progenitors are constructed so that 
they have a constant density inside their half-mass radius. 
Thus one can quickly see from the plot that 'dry' merg- 
ers, with no gas, will produce remnants with densities less 
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than their progenitors and sizes larger than their progeni- 
tors, whereas gas-rich mergers will produce remnants with 
densities higher than their progenitors. The sizes of the rem- 
nants of gas-rich mergers can be similar to or even smaller 
than the sizes of their progenitors. For this set of model 
parameters, constant density evolution occurs at roughly a 
gas-to-stellar- mass ratio of 0.25 (light-blue crosses). 

The stellar-mass FJ plot, in the second panel of Fig- 
ure 2(a) demonstrates a similar effect. Progenitors are not 



shown because they have no values for a, however one can 
see that remnants of mergers with each gas fraction follow 
parallel lines. No rotation is introduced between one gas 
fraction series and the next. 

If elliptical galaxies followed an exact virial relation, 
then one would expect that M star oc a 2 R. Thus plotting 
these quantities against each other gives us a 'virial' projec- 
tion of the fundamental plane. Galaxies following the virial 
relation would fall on a line with a slope of one. Observed 
galaxies do not fall on the expected virial relation. This vari- 
ance from virial scaling is the so-called tilt of the fundamen- 



tal plane. In the third panel of Figure 2(a) we plot M sta r 
versus a 2 R. The line plotted has a slope of one and therefore 
follows a virial scaling. The remnants from each gas fraction 
set fall on the same virial scaling. This demonstrates that 
our model produces no tilt in the fundamental plane for 
remnants when the progenitors' properties are scaled with 
the mass. However, while each gas fraction series follows the 
virial scaling, there is a shift between cases with different 
gas fractions such that larger gas fraction progenitors result 
in remnants with larger ratios of stellar mass to dynamical 
mass. 

In the next experiment we let gas fraction vary as a 
function of mass. We fix the gas mass of the largest progen- 
itor (G3) such that the ratio of gas mass to stellar mass is 
0.25. Then we let the baryonic gas fraction vary as a power 
law with th e baryonic ma s s, MK aB /M barV o ns oc M~J. yons (as 
suggested in iDekel fc Coxi (2006)). We use values of 7 equal 
to 0, 0.5, and 1.0. In Figure 2(b) we show the scaling rela- 



tions of the remnants produced by this series of progenitors. 
Each line of remnants is fixed to the largest mass remnant, 
but the slope of the size-mass and FJ relations clearly de- 
pends on the value chosen for 7. Non-zero values of 7 allow 
for significant rotations in all of the plotted projections of the 
FP. The rotation in the size-mass relation is the direction of 
rotation required if one is to create ellipticals from mergers 
of disks that follow the observed relations. Additionally, the 
fundamental plane relation rotates away from the virial rela- 
tion in the same direction as the observed tilt. For reference, 
the slope of the observed tilt is shown with a dashed line. 
The tilt required is slightly overshot by the 7=1 case, sug- 
gesting that a slightly shallower power law slope would re- 
produce the observed tilt. This result i s compatible with th e 
gas fraction power of 0.7 suggested bv lDekel fc Coxl |2006). 
Thus within our model a gas fraction gradient is capable of 
creating a tilt in the fundamental plane. 

Our model relies on the virial relation to calculate sizes 
and velocity dispersion. However, in our model the central 
dark matter fraction is calculated assuming that no dissi- 
pation occurs within the dark matter halo, and this affects 
our calculation of a. The break from virial scaling results 
from this changing central dark matter fraction, with more 
gas rich progenitors producing a larger difference between 



the dissipational baryons and dissipationless dark matter 
resulting in a lower central dark matter fraction. 

For large values of 7 some curvature is also introduced 
into the scaling relations. Since both the disk and elliptical 
scaling relations are approximately power laws, this puts 
a constraint on the strength of the gas fraction variation 
allowed by the model if we want to reproduce the observed 
scaling relations. However, the expected value of 7 = 0.7 
only produces a modest amount of curvature. 



3.2 Merger Model Parameters 

The merger model contains several parameters with uncer- 
tain values because of the uncertainty in the feedback pre- 
scription used by the merger simulations. Because of this 
uncertainty, the values of the model parameters might re- 
quire adjustment. In order to understand the effect of these 
parameters we vary them systematically and examine the 
results of these variations on the remnant scaling relations. 

We begin by varying the star formation efficiency pa- 
rameter, ei : i. The calibrated value of this parameter (from 
the simulations) is 0.55, but for the experiments below we 
use values of 0.1, 0.5, and 1.0. First, we assume a constant 
gas-to-stellar-mass ratio of 0.5 across the idealized galaxy 
series. Increasing em will increase the number of stars that 
form; however, it has little effect on the remnant scaling re- 



lations (Figure 3(a) I. Specifically, none of the relations ex- 
perience significant rotation as a result of the adjustment. 
For all relations an increase in em simply results in a slight 
shift toward higher stellar mass. 

For the second experiment with em, we introduce a 
mass-dependent gas fraction according to the power law re- 
lation used above, using the expected value of 7 = 0.7. For 
this series, adjusting em produces rotations of the scaling 
relations. Specifically, changing em has a greater effect on 
the mergers with a larger gas fraction, resulting in rotations 
in the size-mass and FJ relations and a slight tilting in the 



FP (Figure 3(b) \ 



Another parameter for which feedback could produce 
some uncertainty is C ra d, which sets the importance of the 
radiative energy term. Since the C* ra d parameter is decou- 
pled from the equation that determines the number of new 
stars, adjusting C ra( j results in no difference in final mass. 
However, an increase of the parameter results in a signif- 
icant reduction of size and increase in velocity dispersion 
for the remnants. The same constant gas fraction and mass- 
dependent gas fraction merger series used for the star forma- 



tion parameter series are plotted in Figures 4(a) and 4(b) 
respectively, with C ra( j taking values of 1.0 (black), 3.0 (red) 
and 5.0 (blue). The value of C ra d determined by fitting to 
the merger simulations was 2.75. As with em, rotation in 
the scaling relations is only seen for the series with a mass- 
dependent gas fraction. Here, a significant rotation is cre- 
ated in the size-mass and FJ relations, but adjusting C ra d 
never introduces a tilt in the FP. This is because the model 
is built on the assumption of the virial theorem, and the 
only portion of the model that violates this assumption is 
the formula for calculating the change in central dark matter 
fraction, which is then used to adjust the velocity dispersion. 
Changing C ra d does not affect this portion of the model. 
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(a) (b) 

Figure 3. Scaling relations for merger remnants produced by the model from our series of idealized progenitors with varying star 
formation parameter e%-i, which takes values of 0.1 (black), 0.5 (red), and 1.0 (blue). 0.55 is the calibrated value. Dotted lines connect 
cases with the same value of e\,\. Progenitors in (a) have a constant gas to stellar mass ratio of 0.5, whereas progenitors in (b) have 
a changing gas fraction with mass (7 = 0.7). The top panels show the size-stellar mass relation of the merger remnants (crosses) and 
progenitors (stars), while the middle and lower panels show the stellar mass FJ relation and the fundamental plane for the merger 
remnants, respectively. The solid line plotted in the lowest panels shows the FP slope for galaxies with a virial scaling. Changing ei : i 
affects the slopes of all three scaling relations, but only when gas fraction varies as a power law of stellar mass (b). 



4 SUMMARY OF OBSERVATIONAL RESULTS 



The Sloan Digital Sky Survey (SDSS) (|York et al.ll2000h has 
provided exquisite statistics on galaxy scaling relations in 
the local universe. IShen et al.l (|2003l ) show that for local 
galaxies the size distribution for each type of galaxy at a 
given stellar mass is log-normal. They provide fitting func- 
tions for the medians of the distributions for both early- and 
late-type galaxies. For late-type galaxies the median (R) is 
described by 



, M \ A M 



(9) 



where 7 = 0.1, a = 0.14, (3 = 0.39, and M = 3 98x IO^Mq 
A com parison of this distribution to that of iBarden et al.l 
(2005), also from SDSS, demonstrates a discrepancy, de- 
spite the fact that the samples contain significant overlap. 
S pecifically, a comparison of the stellar-mass size ridge line 
in lSomerville et all (|2008l ) with the figure in Shen et al., cor- 
recting for the conversion between disk scale length and half 
light radius, shows that the Shen et al. distribution is a fac- 
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(a) (b) 

Figure 4. Scaling relations for merger remnants produced by the model from our series of idealized progenitors with a varying radiative 
energy loss parameter C ra( j, which takes values of 1.0 (black), 3.0 (red), and 5.0 (blue). C ra( j = 2.75 is the calibrated value from the 
simulations. Dotted lines connect cases with the same value of C ra< j. Progenitors in (a) have a constant gas to stellar mass ratio of 
0.5, whereas progenitors in (b) have a changing gas fraction with mass (7 = 0.7). The top panels show the size-stellar mass relation 
of the merger remnants (crosses) and progenitors (stars), while the middle and lower panels show the stellar mass FJ relation and the 
fundamental plane for the merger remnants, respectively. The solid line plotted in the lowest panels shows the FP slope for galaxies with 
a virial scaling. Changing C rac i affects the slopes of the size-mass relation and FJ relation for the cases where gas fraction varies as a 
function of stellar mass (b), but it does not effect the FP. 



tor of ~ 1.5 smaller in radius for a given mass. lDutton et al.l 
find a similar offset between their re-analysis of SDSS 
and the results of Shen et al., and argue that this was pri- 
marily due to their use of circular rather than elliptical aper- 
tures. For this work, we scale the fitting function of Sh en et 
al. to match t he no rmalization of iBarden et al.l (|2005t ) and 
iDutton et ail <|2010| ). 



For early-type galaxies the median radius is described 



by 



.R(kpc) 



M 



(10) 



where a — 0.56 and b — 3.47 x 10~5. The scatter in the 
relation is similar for both late- and early-types with the fit 
to late-type galaxies being 



(CTl - CT 2 ) 



1 + (M/M ) 2 ' 

where <ti = 0.47 and 02 = 0.34. 



(11) 
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The fit for the s tellar mass Faber-Jac kson relation in 
SDSS can be found in iGallazzi et all (2006), and is given by 



log(ff„(km s -1 )) = 0.286 log(M»/M ) - 0.895, 



(12) 



where a v is velocity dispersion and M» is the stellar mass. 
The scatter in the relation is 0.071 dex. 

The fundamental plane relation can be represented in 
a number of projections. However, recent work suggests 
that the tilt of the plane results from a systematic change 
in the central dark matter fraction JZaritskv et al.l 120071 : 
Hopk ins et al]|2008h as the result of a varying effect of dis- 
sipation with mass. A number of observational studies have 
examined the relationship between central dynamical and 
stellar mass, 



M dyn oc M S V 



l+a 
ar 3 



(13) 



and h ave determined a val ue of a sa 0.2 (Pahre et al] 
19981; iGerhard et all l200ll ; iPadmanabhan et all 12004 



Gallazzi et al. 



2006). 



There is little published data on the evolution of the 
FJ and FP to high redshift, though observations sug- 
gest that a ny evolution i n the FP is minimal up 
to z ~ 1 (iHolden et al.l l20ld : ISaglia et aTl I2OI0I : 
iHerbert et al.ll20Hl )7 However, a number of studies have 
examined the evolution of the size-mass relations between 
z ~ and up to z ~ 3 llBarden et al.ll2005l : iMcIntosh et all 
120051 ; iTruiillo et alJ^OOaT The most comprehensive of these 
studies, ITruiillo et ajT |2006). combines data from SDSS, 
GEMS (Galaxy Evolution from Morphology and SEDs) and 
FIRES (Faint Infrared Extragalactic Survey) in order to 
quantify evolution between z — and z — 3 for both early- 
and late-type galaxies. They differentiate galaxies by light 
concentration according to Sersic index n and find that for 
both low-n (late-type) and high-n (early-type) galaxies the 
mean size at given mass at z ~ 2.5 is ~ 2 times smaller than 
today. Specifically they find that at given stellar mass the 



sizes of late- type galaxies evolve as (1 + z) 
the sizes of early-types evolve as (1 + z\ 



3 , whereas 



5 COMPARISON WITH OBSERVATIONS 

We now apply our merger model to major mergers of 
disk-dominated galaxies from the S08 and Millennium 
SAMs. These SAMs produce statistical samples of galaxies 
with properties that are closely matched to those observed 
in the universe. Thus they provide an effective means for 
testing the merger model in a cosmological framework. At 
every time step in the SAMs, we use the data from 
all disk-disk major mergers, regardless of the pre- 
vious or subsequent merger histories of the disks. 
We emphasize that we are making no attempt to 
model subsequent evolution or mergers of the ellip- 
tical remnants. Possible implications of this further 
activity are examined in Section 6. Since the modeled 
merger remnants are a function of the disk progenitor prop- 
erties, we begin by examining the distributions of galaxy 
properties in the SAMs and making comparisons to the ob- 
served distributions. 



5.1 Properties of the Progenitors 

In the model, the most important properties of the progeni- 
tors are initial size, mass, and gas fraction. Thus we begin by 
looking at the size-mass distributions of progenitors in each 



SAM. The relations are plotted for S08 in Figure 5(a) and for 
Millennium in Figure |5(b)| For each figure the progenitors 
are separated into one of six redshift bins according to the 
redshift of the merger. Within each bin, the progenitors 
are divided into mass bins with a width of 0.2 in log(M0). 
The local relations for low-n (solid blue) and high-n (dotted 
red) galaxies are shown for comparison. Additionally, the 
observed redshift ev olution of the size of the low-n galaxies 
ITruiillo et al. 2006) is depicted with the blue-dashed line. 
This is calculated using the median progenitor redshift in 
each redshift bin. 

The size-mass relation produced by the progenitors in 
S08 reproduces the observed relation quite nicely, includ- 
ing evolution with redshift. The Millennium progenitors are 
also fairly close to the observed relation, but in the low- 
est redshift bin they are ~ 50% too large on average. This 
gap lessens with increasing redshift. Also of note is that the 
highest mass bin is typically systematically high. It is also 
interesting to note the difference in slope between the ob- 
served size-mass relations for early- and late-type galaxies. 
As noted before, if a merger explanation of the elliptical 
size-mass relation is to be successful, it must explain the ro- 
tation between the two observed relations. We discuss this 
with respect to our model in the next section. 

Observing cold gas wit hin galaxies is ext remely chal- 
lenging. However a study by iKannappanl (|2004h used a pho- 
tometric estimate of gas fraction in order to calculate the 
changing rat io of gas to s tellar mass (G/S) with galaxies 
from SDSS. ICalura et all (|2007l ) use star formation rates 
from a large SDSS sample combined with theoretical mod- 
eling to estimate (G/S) as a function of both stellar mass 
and redshift. Neither study provides fits to the observed rela - 
tion, but it can be seen from Fig ure 5 of lCalura et ail (|2007l ) 
that, for z < 0.1, 



log(G/5) ~ -0.51og(M/M Q ) + C, 



(14) 



where for blue galaxies C ~ 4.6 and for red galaxies C ~ 
4.1. The values of \og(G/ S) falls approximately in a range 
between -4 and 1.5. 

The relation between G/S and stellar mass in the SAMs 
is shown in Figures 6(a) and 6(b) The observational relation 



(Equation 1141) is also depicted to show that the progenitors 
from the SAMs have gas fraction to stellar mass relations 
with slopes similar to those observed. However, the distri- 
bution of progenitor gas fractions from the SAMs is not ex- 
pected to exactly align with the depicted observed relations 
for two reasons: 1) even the lowest redshift bin (0 < z < 0.3) 
from the SAMs includes redshifts significantly higher than 
those from the observations (z < 0.1), and 2) the progenitors 
are a small subset of all galaxies. Specifically, the galaxies 
shown from the SAMs are those that undergo major merg- 
ers during the specified redshift and are disk-dominated. The 
properties of this subset may differ from both the red and 
blue populations in the observations. In fact, the models of 
S08 have been tuned to reproduce the ga s fraction vs. stellar 
mass observations of Kannappan (2004) for spirals. 

It can be seen that the distributions of G/S in both 
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(a) (b) 

Figure 5. Size Mass relations for the progenitors in the SOS SAM (a) and the Millennium SAM, binned by redshift. Symbols denote 
the median of the progenitor distribution in the SAM, and er ror bars depict the la spread in the distribution. The solid (blue) line and 
the dashed (red) line depict the redshift zero relations from IShen et al.l l|2003h of the low-n and high-n populations res pectively. The 
blue-d ashed and red-dotted lines depict the evolution of the low-n and high-n relations respectively with redshift as given bv lTruiillo et al.l 
(2006). This is calculated using the median redshift for the progenitors in each bin. For this and subsequent scaling relations figures, 
blue and red lines denote observed scaling relations of disk and elliptical galaxies respectively. 



SAMs are bimodal with a gas-rich (blue) and gas-poor (red) 
sequence. While both SAMs capture the slope of the ob- 
served relation, they have gas to baryonic mass power law 
slopes which are significan t ly les s than the 7 = 0.7 value 
suggested in iDekel fe Coxl (|2006T ) and the idealized merg- 
ers presented above. A majority of the mergers are gas-rich; 
99% and 92% of all mergers, for S08 and Millennium re- 
spectively, have total progenitor gas fractions greater than 
0.1. The relation from the SAMs only evolves very mod- 
estly over time, but for both SAMs the fraction of mergers 
that are coming from high gas fraction progenitors increases 
with increasing redshift. This is seen both as a disappear- 
ance of the red sequence (in S08) and a decrease in the typ- 
ical mass of proge nitors as redshift increa ses (both SAMs). 
A recent study (|Tacconi et al.ll2010h found that the 
gas fraction of star-forming galaxies increases by a 
factor of a few as the redshift increases from z — 
to z ~ 2.3. While we have made no selection cuts to 
ensure that our spiral progenitors are star-forming, 
these results suggest that the gas fractions of major 
mergers might increase over time. 

There are also interesting correlations between gas frac- 
tion and size. To demonstrate this, we replot the size-mass 
relation with points drawn for progenitors colored by G/S 



(see Figures 7(a) and 7(b) I. Red points are gas-poor with 
— 2.0 < \og(G/S) < —1.0. Green points have intermediate 
gas fractions with —1.0 < \og(G/S) < 0.0. Blue points are 
gas rich with 0.0 < log(G/S) < 1.0. For both SAMs, at 
given mass there is a significant trend of gas fraction with 
size: the larger radius progenitors have higher gas fractions. 
This is because these galaxies have lower densities and there- 
fore have lower star formation rates and have consumed less 
of their gas, since both SAMs assume a Kennicutt-Schmidt- 
type relationship for star formation (star formation rate den- 
sity is proportional to a power of the gas density). A sim- 
ilar tren d has been observed in a study of nearby 
galaxies (|Catinella et al.ll20ldh . This has important im- 
plications for the evolution of scaling relations via merg- 
ing. Remember that dry merging moves galaxies up and to 
the right on the plot, that is, it increases both mass and 



size (Figure 2(a) I. Because of their location far above the 



late-type scaling relation, dry merging would not take the 
larger radius progenitors on to the desired relation for rem- 
nants. However, since these galaxies are gas-rich they can 
travel with a horizontal or even somewhat downward vector 
when they merge. Additionally, this convergence of merg- 
ing vectors within a single mass bin, such that the gas-poor 
compact galaxies move to larger radii while the gas-rich low 
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(a) (b) 

Figure 6. The relation between G/S and stellar mass for the progenitors in the S08 SAM (a) and the Millennium SAM (b). Points 
are galaxies from t he SAM, The blue (solid) and red (dashed) lines are approximate fits to observations of nearby galaxies in the SDSS 
l|Kannappanll2004l : ICalura et aDl2007h . with blue representing the relation for blue galaxies and red for red galaxies. log(G/5) is plotted 
using a decreasing scale as in the observational studies. The observed relations and SAM progenitors display similar slopes, but the 
relations are not expected to exactly coincide due to differences in redshift ranges and selection effects. 



density galaxies move to smaller radii, results in a reduc- 
tion of scatter for the remnant scaling relations. This effect 
may contribute to the small scatter in the observed early-type 
size-mass relation. 

5.2 Properties of the Merger Remnants 

Given progenitor properties, the merger model allows us to 
predict the sizes, stellar masses, and velocity dispersions of 
the merger remnants. We now exploit this ability to explore 
the buildup and evolution of the remnant scaling relations 
as predicted by the models. In this subsection we use the 
simplified merger model described in Section 2. For compar- 
ison we will show results on the size-mass relation using the 
more elaborate C08 model and summing over merger orbits 
in subsection 5.3 below. 

The remnant size-mass relation for the SAMs is shown 
in Figures [8(a)1 and [8(b)] Note that we are comparing the 
size-mass relations for all observed early-type galaxies at a 
given epoch with that for just the remnants of disk ma- 
jor mergers that occurred in that redshift bin. This ne- 
glects early-types that formed at earlier epochs, and have 
not had a recent merger, as well as early-types formed by 
mixed-morphology mergers. In addition, we only model bi- 
nary mergers, and do not account for multiple mergers or 



subsequent evolution following the formation event. 

Therefore, this may not be a fair comparison, but it is the 
best we can do until we implement the merger model within 
SAMs. 

The model applied to the S08 SAM overpredicts the 
typical observed sizes for a given mass by about 0.3 dex at 
low masses, and the slope is steeper than observed. However, 
the S08 results roughly capture the magnitude of observed 
evolution in average sizes from smaller s izes at high redshift 
to larger sizes at low redshift (see als o iTruiillo et alj|2006l ; 
ITruiillo et ai]|2007l ; iFranx et al.ll2008T ). The slope of the size- 
mass relation is slightly shallower than the observed slope 
at all redshifts, but the merger model has produced a signif- 
icant steepening from the progenitor size-mass relation as a 
result of the mass dependent gas fraction in the progenitors. 

The average predicted size for early-type galaxies in the 
Millennium SAM is too large, sometimes by more than 0.5 
dex. Millennium produces a different normalization for the 
size-mass relation than S08 for two reasons: 1) the progen- 
itors are larger, and 2) the progenitors have less gas. The 
slope of the relation is also much shallower than observed. 
This is likely the result of differences in progenitor gas frac- 
tions between the two SAMs. The evolution in the Millen- 
nium size-mass relationship is such that in the highest red- 
shift bin (2.5 < z < 3.0) the average size at a given mass 
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(a) (b) 

Figure 7. Size-Mass relations for the progenitors in the S08 SAM (a) and Millennium SAM (b), binned by redshift and including points 
for progenitors colored by gas fraction. Red denotes —2.0 < log(G/S) < —1.0. Green denotes —1.0 < log(G/S) < 0.0. Blue denotes 
0.0 < log(G/S) < 1.0. Other symbols and lines are as in Figure [~5(a)| Both SAMs produce systematic gradients in disk galaxy gas 
fractions as a function of size at given mass. This results from higher star formation rates, and therefore faster gas consumption in denser 
galaxies. When dissipative effects are included in a merger model, these systematic gradients result in smaller scatter in the merger 
remnant size-mass relation. 



is ~ 3 times smaller than in the lowest redshift bin. This is 
within the errors on the observed evolution. 

In addition to sizes, the merger model predicts values for 
velocity dispersion (a), so we also examine the stellar-mass 



Faber- Jackson relation (FJ) (see Figures 9(a) and 9(b) I. The 



S08 model roughly produces the correct slope for the local 
FJ relation. However, Millennium produces a slightly steeper 
slope than observed in most redshift bins, particularly at 
high mass. S08 slightly underpredicts (~ 0.05 dex) velocity 
dispersions at z = 0, whereas the offset for the Millennium 
relation is larger (~ 0.2 dex). Since remnants from Millen- 
nium are systematically too large, they also have systemati- 
cally low velocity dispersions. Both SAMs show an evolution 
of the FJ relation with redshift. Within the model there are 
two possible mechanisms for evolution in a at a given mass: 
evolving size and evolving dark matter halo properties. For 
S08 a at a given mass increases by roughly a factor of 1.4 
between z — and 2 = 3; for Millennium a increases by 
a factor of two. These values are commensurate with the 
decrease in size at a given mass for S08 and Millennium, 
suggesting that changing dark matter properties have lit- 
tle effect on the evolution. It is not known yet whether this 
evolution is consistent with the real universe. 

Finally, we examine the fundamental plane (FP) rela- 



tion predicted by the merger model. The model uses the 
virial relation to predict values for a, however the mass used 
in the model is a dynamical mass that includes a contribu- 
tion from dark matter. Therefore, if a tilt is produced by the 
model in the virial projection of the FP, then it is purely 
a result of changing dark matter content in the galaxy cen- 



ters. In Figures 10(a) and 10(b) we plot the relation between 
a 2 r oc Mdyn and stellar mass. The merger model applied to 
the S08 progenitors does not reproduce the observed tilt in 
the FP, with slopes roughly the same as would be expected 
from the virial relation. For Millennium the FP is tilted in 
the same sense as is suggested by observations, with dy- 
namical mass increasing faster than stellar mass, and the 
tilt produced is only slightly less than observed. This tilt re- 
sults from the contrasting dissipational and dissipationless 
evolution of the baryonic and dark matter components com- 
bined with the gas gradient in the progenitors. Lower mass 
progenitors have higher gas fractions and therefore produce 
more compact stellar remnants with smaller central dark 
matter fractions. Higher mass progenitors have smaller gas 
fractions and therefore end up with a sparser stellar center, 
resulting in larger central dark matter fractions. 

Plotting a versus r gives a nearly face-on view of the 
FP, and allows us to determine the portion of the FP be- 
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(a) (b) 

Figure 8. Size-Mass relations for the remnants in the S08 SAM (a) and Millennium SAM (b), binned by redshift. Lines are observed 
relations as in Figure [5(a)| The S08 SAM captures most of the steepening that occurs between the observed late- and early- type scaling 
relations, whereas the slope produced by the Millennium SAM remains near that of the progenitors (disks). Both SAMs show roughly 
the right amount of evolution in size between z=3 and z=0, but both also overpredict the sizes for a given mass. 



ing populated at a given redshift by major mergers of disk 
galaxies. We show this relationship for both S08 (Figure 



11(a)! and Millennium (Figure 11(b) I. Both SAMs show a 



correlation between a and r such that galaxies with larger 
sizes also tend to have larger velocity dispersions. This is 
similar to a correlat ion observed at low redshift in the SDSS 
jGraves et al.ll2009r ). Most interesting perhaps is the notice- 
able evolution with redshift across the face of the FP. At 
high redshift the mergers produce primarily high a galax- 
ies with relatively small sizes. As time progresses toward 
z = 0, the typical a of a merger remnant decreases and the 
typical size increases, as the correlation between the size 
and velocity dispers ion marches across the face of the plane. 
iGraves et all ([2009) also see this correlation between a and 
age, but not between size and age. Also, at any given redshift 
there is a stronger correlation between a and r in our model 
than seen in SDSS. However, as discussed earlier, in real- 
ity each bin would also contain remnants that were formed 
at earlier times, which would weaken the correlation. Fur- 
thermore, one would expect that the remnants from higher 
redshifts would be modified by further merging, which could 
also increase the radius and scatter across the plane. An in- 
teresting question that remains unanswered is whether there 
are specific regions of the FP space that are not populated 
by major merging. We are working on further studies of the 
age and metallicity predicted by our merging model imple- 



mented self-consistently within the SAMs, and will report 
the results in a forthcoming paper (Porter et al. in prep). 



5.3 Scaling relations produced with the C08 
Model 

The new simplified merger model described above captures 
much of the behavior of the C08 model. However, when ap- 
plied to the SAMs, the results of the two models are not iden- 
tical. To illustrate the differences, we show the scaling rela- 
tions produced using the C08 model (Figures ll2U13l and !14[) . 
Aside from using the new, extended model for central dark 
matter fraction (Appendix A), the model is implemented 
as described in C08. For both SAMs the C08 model results 
in a more pronounced steepening of the size-mass relation. 
In fact, at low redshift, for S08 the relation produced has 
a slope that is too steep, whereas for Millennium the slope 
is close to that observed. Of the three scaling relations, the 
size-mass relation is the one that shows the largest differ- 
ence between the two models. However, C08 also produces 
a stronger tilt in the FP relation for both SAMs, with both 
SAMs having a tilt that is close to that observed. Though 
we believe these differences arise from variations in the gas 
fractions and concentrations of the progenitors, we will re- 
frain from a more detailed comparison until the models are 
directly implemented within the SAMs. 
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(a) (b) 

Figure 9. Faber- Jackson relations fo r the remnants in th e S08 SAM (a) and Millennium SAM (b), binned by redshift. The red line is 
the observed relation at low redshift (Gallazzi et al. 2006|). The S08 SAM produces the correct slope for the FJ relation, whereas the 
Millennium SAM relation is somewhat steeper than observed. Both SAMs also underpredict velocity dispersions for a given mass at z=0, 
though the problem is much more severe with the Millennium SAM. Both SAMs show evolution in the normalization of the FJ relation 
since z=3. This is not yet constrained by observations. 



6 DISCUSSION 

A continuing question in the field of galaxy formation is 
the role that mergers play in constructing elliptical galax- 
ies. Sli en et al.l l|2003h have suggested that it is difficult to 
get the proper elliptical size-mass relation from repeated 
major mergers of disk galaxies because of the rotation be- 
tween the early- and late-type relations. Recent SAMs con- 
firm this result, and also produce a scatter i n the remnant 
size-mass relation tha t is much too large (|Shankar et al.l 
l20ld ; lGuo et aUl201Ch . However, these models do not take 
into account the effects of dissipation and a possible de- 
pendence of gas fraction on progenitor mass, as has been 
done in s everal high-resolution hydrodynamical sim- 
ulations dDekel fc Coxll2006l ; I Robertson et al. 2006; 
Hopkins et al.l 120081 ). Here we demonstrate that if the 
effect of dissipation is included within SAMs and if the 
progenitors have a decreasing gas fraction with increasing 
mass, then this results in a steeper size-mass relation for 
the merger remnants compared to the disk progenitors. 
Furthermore, since less dense disks have lower star forma- 
tion rates, the disk galaxy progenitors have increasing gas 
fractions with increasing size at a given mass. Applying a 
dissipative merger model to progenitors with this trend in 
gas fraction also reduces the scatter in the size-mass relation, 



as explained at the end of Section 5.1. In fact, our model, 
with a typical dispersion of 0.2 dex in the size-mass rela- 
tion, has les s than the . 3-0.4 dex dispersion in the relation 
observed by IShen et all l|2003h . That our dispersion is too 
small might be expected, since we are only considering ma- 
jor mergers of disk galaxies, a subset of the total early-type 
population. As noted above, we are also separating 
the size-mass relation according to the redshift of 
the merger. In order to directly compare to obser- 
vations, we would need to include all the mergers 
that have occurred at any given redshift. Comparison 
of the progenitor and remnant properties suggests that the 
overall normalization of the early-type size-mass relation is 
determined by a combination of the size-mass relation of 
the progenitor disk galaxies and the normalization of the 
gas fraction relation. 

In order to capture the observed Faber- Jackson (FJ) 
relation, it is necessary to track the changing central dark 
matter fraction. The results from our model suggest that 
the FJ relation may evolve modestly over time such that 
early-type galaxies at a given mass have larger velocity dis- 
persions at higher redshift. The merger model applied to 
progenitors from the S08 SAM results in a slope close to the 
virial slope, which is slightly shallower than the observed 
slope. The model applied to progenitors from the Millen- 
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(a) (b) 

Figure 10. Fundamental plane plotted as M s t ar versus M^yn for the remnants in the S08 SAM (a) and Millennium SAM (b), binned by 
redshift. The solid red line shows the observed scaling of M^ yn oc M^ r , and the dotted red line shows the virial scaling. The black line 
is a fit to the SAM remnants with M^ yn <x Afi^* and 1 + a is shown on the figure. Neither SAM completely captures the observed tilt 
of the FP. S08 produces tilts that are somewhat less than virial, while Millennium is much closer to the observed tilt with 1 + ot = 1.12 
at z=0. 



nium SAM reproduces the observed tilt in the fundamental 
plane (FP), but this agreement appears to be a fortuitous 
cancellation between sizes that are too large and velocity 
dispersions that are too small at a given mass. The tilt pre- 
dicted by the model is the result of a changing central dark 
matter fraction with mass that results from a mass depen- 
dence of gas fractions in the progenitors. Higher mass pro- 
genitors have less gas and therefore less dissipation during 
their mergers, resulting in remnants with a less concentrated 
stellar center and consequently a higher dark matter frac- 
tion. In the face-on view of the FP, the remnants evolve from 
higher a and lower r at high redshift to lower a and higher 
r at low redshift. 

It is important to note that, since we calculate rem- 
nant properties during post-processing, the scaling relations 
derived using the model in this paper do not capture the 
entire elliptical population at each redshift. Rather they de- 
pict the scaling relations of elliptical galaxies being added 
to the population via major mergers of disks at a given 
redshift. One would expect that the high-redshift galax- 
ies would undergo further accretion and merging, including 
minor merging and dry merging in some ca ses (see e.g. 
iBezanson et al.l [2009 ; lOser et all l2010h . In partic- 
ular, minor mergers between a large, gas-poor el- 
liptical and several small spirals have been shown 



to greatly increase the effective radius of the el- 
liptical whil e leaving the velo city dispersion nearly 
unchanged (Naa b~et al.l [20 09"). Moreover some early- 
type galaxies may be formed via mergers involving one or 
more early-type progenitors, which are also neglected here. 
We find that approximately one-third of the low-redshift 
spheroids at all masses are accounted for by disk-disk major 
mergers, so it could be misleading to compare our model 
outputs with all observed early-type galaxies at various red- 
shifts. 

We also note that there are differences between the re- 
sults of the simple model used here and the previous model 
of C08. Without fully implementing the models within the 
SAMs and tuning the various model parameters, it is dif- 
ficult to know which model will ultimately produce more 
favorable results. The simple model requires fewer unknown 
parameters, and therefore is preferred if it produces rea- 
sonable results. However, the key conclusions of this work 
are that the steepening and reduced scatter of the size- 
mass relation from that of disk galaxies, and the tilt 
in the fundamental plane can all be produced by accounting 
for dissipation during galaxy mergers. While the strength of 
these effects will vary by model, the effects themselves are 
expected to be general features of any realistic merger model 
that accounts for dissipation. 
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(a) (b) 

Figure 11. Nearly face-on view of the fundamental plane plotted as r versus a for the remnants in the S08 SAM (a) and Millennium 
SAM (b), binned by redshift. Both SAMS show some positive correlation between velocity dispersion and size at all redshifts. From z=3 
to z=0 remnants also gradually migrate from higher dispersion and smaller size to lower dispersion and larger size. 



In future work, it will be useful to a incorporate a dis- 
sipational merger model within semi-analytic models. This 
will allow a more comprehensive exploration of the evolu- 
tion of galaxy scaling relations via merging by allowing the 
inclusion of a wider variety of progenitor types, and enable 
the study of correlations with other quantities such as age 
and metallicity. However, even the simple external process- 
ing of SAM progenitors as carried out in this study shows 
that mergers are a viable mechanism for the production of 
a large fraction of elliptical galaxies. Furthermore, we find 
that the elliptical scaling relations are plausibly explained 
given the combination of progenitor scaling relations and 
gas fractions. 
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(a) (b) 

Figure 14. The fundamental plane produced by the C08 model combined with the new model for central dark matter fractions, applied 
to S08 (a) and Millennium (b). The tilt in the fundamental plane is larger than created by the new simple model and for both SAMs it 
is comparable to the observed tilt. 
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7 APPENDIX A: NEW MODEL FOR 

CENTRAL DARK MATTER FRACTION 

In C08 we presented a formula for calculating the dark mat- 
ter fraction in the central part of the merger remnants. We 
define the dark matter fraction with a given radius as 



/d 



M d 



(M dm + M stars ) ' 



(15) 



where M dm and M star s are the dark matter and stellar 
masses inside that radius, respectively. For the merger sim- 
ulations we found that the dark matter fraction inside half 
of the stellar half mass radius could be well-approximated 
by the following formula: 



fdr. 



M d 



+ M d 



M dm ,l + M dm , 2 + C s tars(Ml + M 2 + Mnew ) ' 



(16) 



M dm ,i and Md mi 2 are the dark matter masses inside half of 
the three-dimensional stellar half-mass radii of the progen- 
itors, Mi and M2 are the total stellar masses of the pro- 
genitors, and M ncw is the total mass of stars formed during 
the merger. This expression simply assumes that the inner 
region of the remnant contains the same amount of dark 
matter as the sum of the inner regions of the progenitors, 
and that a fixed fraction, C st ars, of the final stellar mass 
is inside one-half of the three-dimensional stellar half-mass 
radius. 

While this formula works well for the simulated rem- 
nants, it fails outside the regime where Rprogenitor ~ 
Rremnant- One of the benefits of the formula is that since 
it puts all of the central dark matter from the progenitors 
into the central portion of the remnant it allows contraction 
of the halo as a result of the baryonic dissipation. However, 
this contraction is precisely what causes the problem when 
the initial and final radii differ greatly. Specifically, for high 
gas fraction progenitors with large radii, the remnants can 
have much smaller radii and the previous formula predicts 
extreme contraction of the dark matter. 

In order to increase the range of applicability of the 
merger model, we introduce a simpler, more physically intu- 
itive model for predicting the central dark matter fraction 
that does not suffer from the deficiencies mentioned above. 
We begin with the dark matter halo masses and half-mass 



radii, and assume that the two halos merge dissipationlessly 
so that 



(Mi, dm + M 2 , dm ) 2 



It 



•dm,f 



Mj dn 



+ 



M 2 2 dn 



(17) 



This equation can be solved for the final halo half-mass ra- 
dius (-Rdm.f). We use the final radius and mass to fit an 
isothermal profile to the final halo. Then we calculate the 
mass expected inside the stellar half-mass radius and this 
value is used to calculate the final central dark matter frac- 
tion. For this paper we calculate the dark matter fraction 
inside the stellar half-mass radius and use this to compute 
the velocity dispersion of the merger remnant. In addition 
to extending the range of validity, this new method removes 
one of the least certain parameters from the merger model. 
A check against the merger simulations demonstrates that 
the new method produces a larger scatter between predicted 
and measured velocity dispersion, with the fractional rms 
scatter increasing from 0.24 to 0.35. However, most of the 
additional scatter results from variations in orbit, and in this 
work we find that the distribution of orbits found in N-body 
simulations is such that orbital variation plays a minor role. 



